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Spectral function of few electrons in quantum wires and carbon nanotubes as a 

signature of Wigner localization 

Andrea Secchi^'-^ and Massimo Rontani^'El 

^CNR-NANO S3, Via Campi 213a, 41125 Modena, Italy 
^ Dept of Physics, University of Modena and Reggio Emilia, ttaly 
(Dated: March 20, 2012) 

We demonstrate that the profile of the space-resolved spectral function at finite temperature 
provides a signature of Wigner localization for electrons in quantum wires and semiconducting 
carbon nanotubes. Our numerical evidence is based on the exact diagonalization of the microscopic 
Hamiltonian of few particles interacting in gate-defined quantum dots. The minimal temperature 
required to suppress residual exchange effects in the spectral function image of (nanotubes) quantum 
wires lies in the (sub-) Kelvin range. 
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After half a century of research, electrons in one di- 
mension still attract attention as a paradigm of inter- 
acting behavior that deviates from Fermi liquid the- 
ory, exhibiting e.g. spin-charge separation or solid-like 
order ji"— This impulse comes from recent experiments 
in systems with high aspect ratio — cleaved-edge over- 
growth structures f-"-- multiple-gate quantum wires,^ car- 



tively one-dimensional as their transverse and longitudi- 
nal degrees of freedom are decoupled. The refinement 
of such devices allows to easily reach the dilute regime 
of electron density yet minimizing the impact of disor- 
der. At sufficiently low density, the Coulomb energy gain 
overcomes the kinetic energy cost of localization, hence 
electrons are expected to freeze their motion in space 
forming a regular array — a Wigner correlated solidj ^^'^^ 

So far, the evidence of Wigner localization has ulti- 
mately relied on the measure of the energy gap between 
ground and low-lying excited states, which vanishes in 
the dilute limit i^i^ '^^i^^ This excitation energy decreases 
gradually from the liquid- to the solid-like regime, as an 
effect of both quantum fluctuations and samples' finite 
size — systems often act as quantum dots (QDs) in the 
Coulomb blockade regimei^i^""— Therefore, an alterna- 
tive signature of the electron solid, directly related to the 
wave function, would be desirable. A possible observable 
is the momentum-resolved spectral function (SF) — the 
quasiparticle wave function square modulus in recipro- 
cal space.— Unpromisingly, it was predicted that the 
SF was qualitatively similar in both Wigner and non- 
interacting limitsi^ and that any distinctive structure of 
the SF was washed out by temperature^^ 

In this Communication we demonstrate that the space- 
resolved spectral function of few electrons provides a clear 
signature of Wigner localization at temperatures above 
Tex, that is the characteristic scale of exchange inter- 
actions. This fundamental observable may be accessed 
through scanning tunneling spectroscopy (STS)i^i2i Our 
exact diagonalization?^"^^ (ED) results show that the SF 
resembles the charge density for T > Tex, displaying 
TV peaks as the TVth electron tunnels into a Coulomb 
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FIG. 1. (Color online) (a-c) SF signal 5 vs x at various T 
for TV = 2 ^ iV = 3. (a) Quantum wire (QW). The length 
unit is Iqo ~ 23.9 nm. (b) Carbon nanotube no. 1 (CNTl). 
^QD = 26.8 nm. (c) Carbon nanotube no. 2 (CNT2). ^qd = 
17.9 nm. (d) SF signal S vs 1 at T = 10 K for a CNT as 
hojo is varied. The five curves, from red [light gray] to black, 
correspond to fiwo = 5, 20, 60, 100, oo meV, respectively. The 
length unit is ^qd and the remaining parameters are those of 
CNT2. (e) Tex vs &jo for QW, CNTl, and CNT2. The filled 
circles are the values pertinent to panels (a-c). 



blockaded QD already containing iV — I electrons. The 
peak-to- valley ratio of such image allows to assess directly 
the onset of Wigner localization. In sharp contrast, for 
T <C Tox the SF is system-dependent and unrelated to 
N. Overall, the joint measurements of N and of the SF 
are able to unveil the Wigner solid. 
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The scenario reported in this Communication agrees 
with the theory of the 'spin-incoherent' Luttinger 
hquid.— On the other hand, crucial approximations of 
this theory poorly reproduce key experimental features of 
systems with a moderate number of electrons and 
most noticeably: (i) finite-size effects are prominent}^"— 
(ii) the occurrence of the band gap^ and band 
curvaturSi^!^ may not be neglected. 

In our ED approach we take into account all many- 
body correlations as well as the effects of spin-orbit cou- 
pling, valley degeneracy, band curvature through the ef- 
fective mass We assume the QD confinement po- 
tential along X to be harmonic, Vqj:,{x) = moj^x^ /2^^^ 
since this is the generic low-energy form for gated QDs 
embedded in quantum wires (QWs) and semiconducting 
carbon nanotubes (CNTs), as in Refs. [sl. I8l-[l2i We con- 
sider N electrons in a QW interacting through a screened 
Coulomb interaction l/(a;, x') = e^e~^[(a; — x')^-|-A^]~^/^, 
with A being a short-range cutoff and e the dielectric con- 
stant. The CNT Hamiltonian is more complex, due to the 
presence of valleys K and K', spin-orbit coupling, inter- 
and intra-valley interactions. The interaction potential 
interpolates between Coulomb and Hubbard-like behav- 
ior: see Ref. [2^ for details. In all cases we diagonalize 
the Hamiltonian in the space spanned by the Slater de- 
terminants built by filling with N electrons the lowest 60 
spin-orbitals in all possible ways. 

The outcome of the ED consists in the ith excited iV- 
body states, |-/V, and their energies . The SF for a 
given initial state at T = 0, Ais!-i^i{x,u}), is 
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FIG. 2. (Color online) QW q{x) of A'' electrons (black curve) 
and SF at T = for the — 1 — !> A'^ ground-state transition 
(solid red [gray] curve) vs x for A'^ = 2 (a), 3 (b), 4 (c), 5 
(d). The dashed curves are the non-interacting SFs, rescaled 
by a factor 1/6. AN-ifi{x,uj) is integrated over a tiny energy 
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The length unit is i?QD = 23.9 nm. 
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with ^''{x) being the operator creating an elec- 
tron at X and fiw its resonant tunneling energy, 
whereas the ground-state charge density is g{x) — 

N-^ {N, 0| *t(a;)^(a;) \N, 0). Ideally, the STS differential 
conductance at vanishing temperature and bias, dl/dV, 
is proportional to AN^i^Q^x^SF/fi), with sp being the 
Fermi energy of the STS tip.^'^ If the thermal broadening 
ksT is larger than typical energy spacings,'^'^ dl/dV is 
proportional to 



g — duj 



duj 



At{x,uj), 



(1) 



where /(•) is the Fermi distribution function and 
Ax{x,uj) is the finite-temperature SF, 

At{x, ^) = ^ E e^^^'""^iv-i,.(x, u), (2) 



with 13 = l/iksT) and Z = cxp (-/S-Sf -^). In the 
following we tune ep appearing in /(•) to match the A^ — 
1 ^ A^ transition between ground states. 

To immediately grasp the key results of this Commu- 
nication, consider in Fig. [IJc) the SF signal Q induced 



by the tunneling of the third electron into a realistic 
CNT (here labeled CNT2) at various temperatures. At 
low temperature, T <^ Tex — 0.017 K, the signal shows 
two peaks that originate from the symmetries of the 
ground states. Above Tex the signal shows three well 
resolved peaks that resemble the charge distribution of 
three Wigner- localized electrons [cf. g{x) in Fig. IH^b)]. 
The peak-to-valley ratio at T 3> Tex directly measures 
the degree of spatial localization. This is apparent in 
Fig. md), as the impact of few-body correlations is re- 
duced by increasing the QD confinement energy Hujq and 
hence the ratio of kinetic to Coulomb energy. Whereas 
for strong correlations (red [light gray] curve) the three 
peaks are well resolved since electrons separately localize 
in space, as the interaction is turned off Q becomes fea- 
tureless (black curve for HijJq = oo), clearly discriminating 
between Wigner and weakly-interacting regimes. 

Also the magnitude of Tex points to electron corre- 
lation, the lower the temperature Tex, the stronger the 
localization. Tex varies significantly for typical QWs and 
CNTs as a function of device parameters, like huiQ. In 
Fig. [TJd) the increase of hujo quenches correlations and 
amplifies the effects of Fermi statistics, raising Tex. For 
realistic parameters [cf . circles in Fig. [TJd)] , we find that 
Tex may be as low as 10 mK in some semiconducting 
CNTs. On the other hand, Tex is one-two orders of mag- 
nitude higher in quantum wires (QWs), Tex 1 K, as an 
effect of the different impact of screening. 

Figure [5] shows g{x) and the zero-temperature SF up to 
five electrons for a typical wire QD. In the ED we chose 
A = 5 nm, bulk GaAs parameters, and huo = 2 meV as 
single-particle energy spacing, providing a characteristic 
QD length ^qd = (fi./mwo)'^^^ = 23.9 nm. As it is seen 
from the spread of g{x) along the axis (black curves in 
Fig. HI), this corresponds to a typical size of ~ 200 nm for 
A^ = 5, which is comparable to the size Lioc — 230 nm of 
the QD in the dilute limit of Ref. H (Table I). The profiles 
of g{x) point to the partial localization of the A^ electrons 
as N peaks emerge from a featureless liquid droplet. On 
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FIG. 3. (Color online) CNT g{x) of electrons (black curve) 
and SF at T = for the N — 1 N ground-state transition 
(solid red [gray] curve) vs x for A'^ = 2 (a), 3 (b), 4 (c), 5 
(d). The dashed curves are the non- interacting SFs, rescaled 
by a factor 1/6. An-i,o{x,ijj) is integrated over a tiny energy 
range around — Eq~^ . The length unit is ^qd = 26.8 nm. 



the other hand, the SFs for the N — 1 ^ N ground state 
transitions (solid red [gray] curves in Fig. [2]) are qualita- 
tively distinct from g{x) at given N. This is patent for 
N = 3,4, 5, with SFs displaying one, one, and two nodes, 
respectively, whereas the corresponding £'(a;)'s have no 
nodes. As seen in Fig. [IJa), Q is insensitive to tempera- 
ture in the range T <C Tex =1.1 K, with fc^Tex being the 
energy splitting between the lowest two-electron triplet 
and singlet states. Here fcsTcx provides the energy scale 
of low-lying excitations. 

The SFs appearing in Fig. [2] are similar to those ob- 
tained in the absence of interaction (dashed lines). In 
the non-interacting limit the SF is the square modulus 
of the orbital occupied by the A^th electron that enters 
the QD.~ Such orbital has zero, one, one, and two nodes 
for iV = 2, 3, 4, 5, respectively, since electrons fill in each 
orbital level twice due to Kramers degeneracy. As the 
symmetries of the quantum states do not change in the 
considered range of interaction, no qualitative differences 
are seen for the interacting SFsJ^ 

In the following we discuss two examplar CNT cases. 
Figure [3] is the analogue of Fig. [5] for the CNT device 
investigated in Ref. d (here labeled as no. 1). The ED 
parameters (Swq = 8 meV, e = 3.5, radius R — 3.6 nm) 
were chosen in order to reproduce the measured chemi- 
cal potentials, as detailed in Ref. [2^ Apart from length 
renormalization (^qd = 26.8 nm), charge densities (black 
curves in Fig. [3]) are similar to those of the QW (Fig. [5]). 
On the contrary, the zero-temperature SFs for the CNT 
(solid red [gray] curves in Fig. [3]) are drastically differ- 
ent from those of Fig. [2l being all nodeless except for 
the N — 4: ^ N — 5 transition. This trend is qual- 
itatively similar to the non-interacting filling sequence 
(dashed lines in Fig. [3]) , as each CNT level is four- fold 
degenerate in the absence of spin-orbit coupling, due to 
both spin and valley degeneracies. The spin-orbit inter- 
action splits the multiplet into two doublets (here sepa- 
rated by AEso = 0.367 meV) but leaves the spin-orbitals 
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FIG. 4. (Color online) SF signal S vs a; at T = 0.1 K (dashed 
curve) and T = 0.5 K (solid red [gray] and green [light gray] 
curve) for the CNT QD no. 2. The green [light gray] curve is 
the non-interacting case. The following N — 1 ^ N transitions 
are considered: (a) = 2, (b) = 3, (c) A^ = 4, (d) A^ = 5. 
q{x) for the A'^-body ground state (black curve) is plotted for 
comparison. The length unit is ^qd = 17.9 nm. 



unchanged. Since AEso remains the energy scale of the 
low-lying excitations even in the presence of interactions 
in the sample no. 1 (Refs. i and [H), the SFs of Fig. [3] 

are unaffected by temperature for T ^ AEso/ks ~ 4 K 

[cf. Fig.inb)]. 

The SFs of Figs. [Sja), (b), and (c) display A^ peaks as 
the £i(x)'s, and placed approximately in the same loca- 
tions. This genuine effect of interaction, reminescent of 
the partial Wigner localization occuring in the QD, takes 
place also in the QW [Fig.[2l^a)] and it has been observed 
for the tunneling of the second electron into elongated 
self- assembled InAs QDs."^^ 

The second CNT QD that we study experiences 
stronger interactions than the first one, as an effect of 
the smaller energy spacing {huo — 5 meV), dielectric 
screening (e — 2.5), radius {R = 1 nm). In terms of pa- 
rameters, the CNT QD no. 2 lies in the middle between 
the devices investigated in Refs. d and i (see Ref. IH for 
their placement in a phase diagram). As it is shown in 
Fig. U (black curves), the peak-to- valley ratios of charge 
densities are about twice as large as those in Figs. [5] and 
[3l hence electrons are strongly localized. 

In this Wigner regime the low-lying excited states are 
easily thermally populated and are highly degenerate. 
The reason is that exchange interactions between local- 
ized electrons are suppressed (Tex is only 17 mK), there- 
fore each electron may flip its spin independently from 
the others at low energy costi^^ This is true also for the 
isospin — the orbital angular momentum along x labelling 
valleys K and K'. The overall result is that the many ex- 
cited states which differ only in the (iso)spin value are al- 
most degenerate. This situation is illustrated in Fig.lSJa), 
where we depict the excited-state ladders for N — 2 and 
A^ = 3. In the ultimate Wigner limit one expects the 
ground state to be exactly 4^ degenerate, as each elec- 
tron may flip its spin and isospin in four different ways. 
The effect of spin-orbit interaction (here AEso = 1-32 
meV) is to split the ground-state multiplet into equally 
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FIG. 5. (Color online) SFs of the CNT QD no. 2. (a) Lowest 
excited states for N = 2 and N = 3. The numbers label 
sub-multiplet degeneracies, (b) Blow-up of panel a (not in 
scale). Arrows point to allowed tunneling transitions, (c) SFs 
of allowed transitions vs a; at T = 0. The various AN-i,oix, uj) 
are integrated over tiny energy ranges around the resonance 
energies highlighted by the arrows in panel b, according to 
the color [gray tone] code. The length unit is ^qd ~ 17.9 nm. 



spaced sub-multiplets (three and four for = 2 and 
= 3, respectively). 

Residual exchange interactions further split the lev- 
els, as illustrated by the blow-up of the two lowest sub- 
multiplets for iV = 2 and = 3 shown in Fig. [5{b) (not 
in scale). For example, states A and B for = 2, having 
the orbital part of the wave function even and odd with 
respect to reflection symmetry, respectively, by definition 
are split by fcsTex = 1.38 /xeV. Similarly, states a, b, c for 
A'^ = 3 have distinct orbital symmetries and are split by 
4.45 /xeV (a from 6) and 2.26 /ieV (6 from c), respectively. 

It is clear that already at Tex = 17 mK the excited 
states of Fig. [5jb) are significantly populated. Therefore, 
even at low temperatures — say T = 100 mK for state-of- 
the-art STS — the SF is a statistical mixture of several 
transitions, explicited in Fig. [S] The arrows depicted 
in Fig. [Sfb) point to the low-energy transitions between 
N = 2 and A^ = 3 that are allowed by spin, isospin, 
and orbital symmetries.-^ The corresponding SFs at T = 
are shown in Fig. EJc). Each plot of Fig. [5jc) maps 
onto a different transition [the arrow of like color (gray 
tone) in Fig. O^b)] except the red [light gray] plot which 
is identical for both A ^ a and B ^ b resonances. 

The SFs plotted in Fig. [Sfc) differ among themselves 
with regards to both the number of peaks and their rela- 
tive intensities. The variations are dictated by the sym- 



metries of initial and final states involved in the tunnel- 
ing transition. For example, the fundamental transition 
A ^ a between two- and three-electron ground states dis- 
plays two peaks located at opposite positions (red [light 
gray] curve), whereas all other depicted SFs have three 
peaks each. This shows that the SF at T = may greatly 
deviate from the charge-density profile, corroborating the 
findings of Figs. [5] and [31 On the other hand, the loca- 
tions of the maxima of curves in Fig. [Sjc) coincide. 

By statistically averaging the SFs of low-lying transi- 
tions, as those shown in Fig. [5jc), one obtains the finite- 
temperature signal G [cf. Eq. ((!])]. Figure |4] shows the 
pattern of ^ at T = 0.1 K (dashed curves) and T = 0.5 
K (solid red [gray] curves) for transitions A^ — 1 — > A^ up 
to A^ = 5. The small dependence of Q on the temperature 
exhibited in Figs. HJ^c-d) [see differences between dashed 
and solid red (gray) curves] is a finite-size effect due to 
the form of the potential Vqd {x}, as the electron density 
slightly increases with A^ (Ref. l23l ) . Remarkably, Q has a 
regular behavior as a function of A^ already at T = 0.5 K, 
systematically displaying A^ peaks of comparable heights 
whose positions are close to (but non coinciding with) 
the locations of the maxima of g{x) (black curves). 

The pattern of Q shown in Fig. |4] — exhibiting high 
peak-to-valley ratio — is peculiar of the Wigner regime 
and should be contrasted with the featureless, non- 
interacting profile (green [light gray] curves). Indeed, 
the low- lying excited states, as those of Fig. Elja), have 
all roughly the same orbital wave function modulus, sim- 
ilar to the vibrational wave function of nuclei of poly- 
atomic moleculesi^ The differences among orbital states, 
as well as those among SFs [cf. Fig. [SJc)], originate from 
the different nodal surfaces. Since the weight is mainly 
localized around the equilibrium positions of electrons, 
nodeless interstitial regions may hardly be distinguished 
from nodal regions. Therefore, the statistical average of 
excited states shows a regular trend, linked to the posi- 
tions of localized electrons. 

At sufficiently high temperatures, the SF signals Q of 
all investigated devices behave similarly. This is shown 
in Fig. [T] for the tunneling transition A^ = 2 A^ = 3. 
Above their respective temperatures Tex, highlighted as 
circles in Fig. [TJd), all Q profiles exhibit three peaks of 
similar height: see the curves for T = 4, 10, 0.1 K in 
Figs. [TJa), (b), (c), respectively. The different peak-to- 
valley ratios of these curves measure the degree of Wigner 
localization, the lower Tex, the higher the ratio. 

At T Tox the central peak of the QW signal is 
depleted again [red (light gray) curve for T = 30 K in 
Fig. [TJa)], whereas CNT profiles [red (light gray) curves 
in Figs, [lib) and (c)] remain stable well above 50 K. This 
change is due to the excitation of the energy scale associ- 
ated to charge. In fact, the charging energy of the QW, 
estimated as the energy difference between the resonance 
energies of the first two electrons, is 5.1 meV, that is com- 
parable with 30 K, whereas CNT charging energies are 
much larger (~ 20 meV). 

In conclusion, we have shown that the spatial depen- 
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dence of the spectral function provides a clear fingerprint 
for Wigner localization, as the temperature overcomes 
the energy scale of exchange interactions. This temper- 
ature is low enough in both semiconducting carbon nan- 
otubes and quantum wires to make scanning tunneling 
spectroscopy feasible. This effect has not been seen in 
past experiments j^i^ likely due to metallic screening, as 
well as to the presence of disorder and scattering from 
boundaries. We hope our prediction may stimulate fur- 
ther work along this path. 
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